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Abstract. We report on Hubble Space Telescope Faint 
Object Spectrograph observations of the double QSO HS 
1216+5032 AB (z e (A) = 1.455; z e (B) = 1.451, and an- 
gular separation 9 = 9'.'1). The spectral coverage is 910 
A to 1340 A in the QSO rest-frame. An unusual broad- 
absorption-line (BAL) system is observed only in the B 
component: maximum outflow velocity of ~ 5 000 km s ; 
probably a mixture of broad and narrow components. Ob- 
served ions are: H I, C II, C in, N III, N v, O VI, and possibly 
S IV and S VI. We also discuss two outstanding interven- 
ing systems: (1) a complex Civ system at z = 0.72 of 
similar strength in A and B, with a velocity span of 1500 
km s^ 1 along the lines of sight (LOSs; LOS separation: 
S±[z = 0.72] « 75 h^Q kpc); and (2) a possible strong 
Mgn system at z = 0.04 observed in B only, presumably 
arising in a damped Lya system. 

We assume HS 1216+5032 is a binary QSO but dis- 
cuss the possibility of a gravitational lens system. The size 
of Lya forest clouds is constrained using Sj_ « 80 
kpc at redshifts between z = 1.15 and 1.45. Four Lya 
systems not associated with metal lines and producing 
lines with Wo > 0.17 A are observed in both spectra, 
while five appear in only one spectrum. This sample, al- 
though scarce due to the redshift path blocked out by 
the BALs in B, allows us to place upper limits on the 
transverse cloud sizes. Modelling the absorbers as non- 
evolving spheres, a maximum-likelihood analysis yields a 
most probable cloud diameter D = 256 h^Q kpc and 2a 
bounds of 172 < D < 896 h^ kpc. If the clouds are mod- 
elled as filamentary structures, the same analysis yields 
lower transverse dimensions by a factor of two. Indepen- 
dently of the maximum-likelihood approach, the equiva- 
lent width differences provide evidence for coherent struc- 
tures. The suggestion that the size of Lya forest clouds 
increases with decreasing redshift is not confirmed. Fi- 
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nally, we discuss two z a « z e Lya systems observed in 
both QSO spectra. 

Key words: Cosmology: observations - Quasars: individ- 
ual: HS 1216+5032 - Quasars: general - absorption lines 



1. Introduction 

HS 1216+5032 was discovered in the course of the Ham- 



burg Quasar Survey (Hagen et al. 1995 ). It was later con- 
firmed as a double QSO at z e — 1.45 by Hagen ct al. 
( |1996| ) through direct images and low-resolution (FWHM 
« 20 A) spectroscopy in the optical range. The B magni- 
tudes of the bright and faint QSO images (hereafter "A" 
and "B" , respectively) are B A = 17.2 and B B = 19.0. The 
separation angle between A and B is 6 — 9'.'1, and the 
emission redshifts are z e (A) = 1.455 and z e (B) = 1.451. 
From the available optical data of HS 1216+5032 AB there 
appears to be some evidence favoring its physical-pair na- 
ture instead of a gravitational lens origin (this subject will 
be discussed in §3.1), so the projected distance^ between 



LOSs as a function of redshift is known (78 h 50 kpc at 

This paper presents ultraviolet (UV) spectra of QSO A 
and B taken with the Faint Object Spectrograph (FOS) 
on-board the Hubble Space Telescope (HST). The prime 
goal of these observations was originally to investigate the 
Lya forest in the little explored redshift range between z = 
0.8 and z — 1.4. In particular, the projected separation 
between the LOSs to HS 1216+5032 A and B is convenient 
because it samples well the expected Lya cloud sizes of 
hundreds of kpc. 

For the redshift range accessible from the ground, a 
pioneering work on cloud sizes was made by Smette et al. 
( |j~992T ) using the gravitationally lensed double QSO UM 
673 (9 = 2'.'2) with known lens geometry. All observed Lya 
lines were observed to be common to both spectra and a 
statistical lower limit of D > 12 h§ Q kpc could be set 



1 Throughout this paper we assume an (Qq, Ao) = (1,0) cos- 
mology and define hso = Ho/(50 km s~ 1 Mpc~ 1 ). 
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for the transverse sizes (spherical clouds, no evolution). 
A similar result (Smette et al. 1995) was found for HE 
1104-1805 AB (0 = 3") but here the position of the lens 
z\ was unknown. A transverse size of D > 60 kpc was 
estimated for 2|»1, 

Surprisingly, QSO pairs at very large projected dis- 
tances seem to still show Lya lines common to both spec- 
tra. In the spe ctra o f LB 9605 and LB 9612 (0 = 1.'65), 
Dinshaw et al. ( 1998 ) find 5 such lines within 400 km s _1 



and derive a most probable diameter of 1520 h^ kpc at 
1 < z < 1.7. If the lines arise in different absorbers, how- 
ever, an upper limit of ~ 1140 kpc is derived. Fur- 
ther studies in the optical range have been carried out 
by Dinshaw et al. fll994j; 6 = 9"5, D > 160 kpc at 



1.7 < z < 2.1) and Crotts et al. (|l994j ; same results). 

At low redshift (z 55 1.5), Lya clouds show a very flat 
redshift distribution: only 16 ±2 Lya lines with Wq > 0.32 
A are expected in the wavelength interval covered by the 
FOS on the HST (Weymann et al. |1998j ). Therefore, the 
size-estimate uncertainties for low-redshift Lya absorbers 
are large. Studies using FOS spectra of QSO pairs have 



been made by Dinshaw et al. (1997; 9 = lf29, D > 600 
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kpc at 0.5 < z < 0.9), and more recently by Petitjean 



et al. (|1998t 9 = 36", D ~ 500 kpc at 0.8 < z < 1.4) 



All the aforementioned limits on Lya cloud sizes result 
from rather crude models which assume non-evolving and 
uniform-sized absorbers. Such simple models might be un- 
realistic, as suggested for example by a trend of larger size 
estimates with larger LOS separation (Fang et al. 1996) 
or by the claim (Fang et al. 1996, Dinshaw et al. 1998) 



that the cloud size increased with cosmic time. However, 
we note that none of the above claims is confirmed by 
D'Odorico et al. ( 199£ ), using a larger database. Clearly, 
the issues of geometry and size of Lya absorbers remain 
open due partly to a lack of size estimates at low redshift. 

In this paper we use a likelihood analysis (§ ^) that 
attempts to constrain the size of Lya clouds in front of 
HS 1216+5032 A and B based on simple considerations 
on shape and evolution. We model the absorbers either as 
non-evolving spheres or as filaments. The technique takes 
advantage of the information provided by line pairs ob- 
served in both spectra at the same redshift, hereafter re- 
ferred to as "coincidences" , and Lya lines observed only 
in one spectrum, referred to as "anti-coincidences" . A like- 
lihood function is then constructed using an analytic ex- 
pression for the probability of getting the observed num- 
ber of coincidences and anti-coincidences in each of the 
two geometries. 

In general, modelling the geometry of Lya absorbers at 
low redshift on the basis of line counting is made difficult 
because the line samples are intrinsically small, implying a 
poor statistics. In addition, HS 1216+5032 B shows broad 
absorption lines (BALs) caused by several ions with tran- 
sitions in the UV, so the effective redshift path for the 
detection of Lya lines is even shorter. 



Intervening metal absorption lines of interest are pre- 
sented in § [I| The BAL systems are described in § M. The 
conclusions are outlined in § [?]. 

2. Data reduction 

2.1. Observations and wavelength calibration 

UV spectra of HS 1216+5032 A and B were obtained on 
1996 November 6 with the FOS on board the HST. Tar- 
get acquisition and spectroscopy were done using Grating 
G270H with the red detector and the 3."7 x 3."7 aper- 
ture. This configuration yields a spectral resolution of 
FWHM = 2 A and a wavelength coverage from 2222 A 
to 3277 A (Schneider et al. |1993| ). Total inte gration times 
were 1980 and 10400 s for QSO components A and B, 
respectively. The spectrum of image B resulted from the 
variance- weighted addition of four exposures. The signal- 
to-noise ratios at Lya emission are S/N = 40 (A) and 25 
per ~ 0.5 A pixel, falling down to 15 in the blue part of 
the A spectrum, and to near one at the BAL troughs in 
B. 

The flux-calibrated spectra and their associated 1 a er- 
rors are sho wn in Figure [l]. The dotted lines represent the 
continua (§ ) and the tick-marks indicate the position 
of 3cr absorption lines (§ ||). 

Since we are interested in comparing absorption sys- 
tems along both LOSs, we must check for possible small 
misalignments of the zero-point wavelength scale between 
both spectra. Such differences may arise when the tar- 
gets have not been properly centered in the aperture of 
the FOS, and/or when - unlike here - the science and 
wavelength calibration exposures have not been performed 
consecutively. Galactic absorption lines common to both 
spectra were then required to appear at the same wave- 
length, under the assumption that both LOSs cross the 
same cloud (which should be true for the small separa- 
tion between LOSs at z k 0). There are two Galactic 
absorption lines, Fen A2600 and Mgn A2796, common to 
both spectra and apparently not contaminated with inter- 
galactic lines. For these lines AA(Fe n A2600)= —9 km s _1 
and AA(Mgn A2796)= +4 km s _1 are obtained, which lie 
within the fitting procedure errors (~ 10 km s" 1 ) and 
the FOS limiting accuracy (~ 20 km s _1 ). These small 
differences imply that both spectra are well-calibrated; 
consequently, we made no corrections to the zero-point 
wavelength scale. 



2.2. Emission redshifts 

To derive the emission redshift of HS 1216+5032 A the 
Lya emission peak was fitted with a Gaussian in the re- 
gion between 2959 and 3005 A. Excluding the wavelength 
regions with absorption features in the blue wing of the 
line yields z e (A) — 1.4545 ± 0.0001, where the la error 
comes from the wavelength uncertainty. For QSO compo- 
nent B, a similar analysis is extremely difficult because 
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the Lya emission line is severely distorted by the Lya 
and Nv BALs. Fitting instead the Nv emission peak be- 
tween 3016 and 3065 A, and using the mean rest-frame 
wavelength of the doublet A (Nv)= 1240.15 A, yields 
z e (B) = 1.4509 ± 0.0006. In the following we adopt the 
nominal values z e (A)= 1.455 and z e (B)= 1.451. 

The redshifts of A and B are not the same within mea- 
surement errors, with QSO component B being blueshifted 
by ~ 400 km s _1 relative to A. Different emission redshifts 
are expected in the physical pair hypothesis, but the mea- 
sured difference between HS 1216+5032 A and B cannot 
rule out the possibility that the QSO may still be gravita- 
tionally lensed. This is because two different transitions, 
Lya and N v, are being compared and blueshifts of a few 
hundreds km s _1 between high and low-ionization emis- 
sion lines are common (e.g., Hamann et al. 1997; Tytler 



& Fan 1992). This point is further discussed in § 3.1 



Other prominent emission lines in the spectrum of A 
are: Ovi A1033 (blended with Ly/3) and probably Cm 
A977 and Siv A1062. In the spectrum of B, the BAL 
troughs allow clear identification of only Lya and Nv 
A1240, but the "undulating" shape of the observed flux 
for A < 2600 A suggests that S vi A939, C in A977, N in 
A990, Ovi A1033, and Siv A1062 have broad absorption 
troughs on the blue side of the expected QSO rest wave- 
lengths. 

2.3. Continuum fitting 

Due to line blending in the Lya forest, defining a quasar 
continuum shortward of Lya emission is not trivial at our 
resolution. For B, this difficulty is increased by the BAL 
profiles. A power-law shaped continuum / oc v a , although 
an acceptable representation in the optical range, is defini- 
tively not able to fit even the UV continuum spectrum of 
quasar image A, leaving residuals far larger than 3 a at 
apparently featureless spectral regions. 

We decided to fit the UV QSO continuum for each 
spectrum with cubic splines, using the following simple 
and semi-automated fit algorithm. First, cubic splines are 
fitted through a given number of points (xi,yi) thought 
to represent the continuum, typically 50 in each spectrum. 
The point ordinates t/, are then corrected by an amount, 
such that the mean of the differences between the new or- 
dinate and all the flux values within la (defined by the 
first fit) in a 6 A wide spectral window around Xi equals 
zero. The corrected (xi,yi) are used to perform a new fit. 
This procedure is iterated until it converges to a reason- 
able continuum, provided enough featureless regions are 
available in the data. The fit routine is robust in the sense 
that different start values always lead to the same con- 
tinuum, but the spectral windows need to be carefully 
chosen, especially in the B spectrum. The final continuum 
is shown as dotted line in Fig. [l]. Note that we have not 
attempted to model the QSO intrinsic continuum of B, 
which is evidently distorted by the BALs. Division of the 



flux by this continuum results in the normalized spectra 
of HS 1216+5032 A and B. 



3. Absorption line parameters and identifications 

Tables [l] and [| list absorption lines found at the 3cx level 
in HS 1216+5032 A and B, respectively. Lines in the nor- 
malized spectra were fitted with Gaussian profiles without 
constraints in central wavelength, width or amplitude. In 
the case of obvious blends two or more Gaussian compo- 
nents were fitted simultaneously. 

The fit routine attempts to minimize x 2 between model 
and data, considering la flux errors (flat-fielding, back- 
ground subtraction and photon statistics noise); but the 
uncertainties introduced by the continuum placement are 
not included. 

An inherent feature of % 2 minimization is the non- 
uniqueness of the solution due to the eventual presence of 
more than one minimum in the parameter space. In those 
cases, the fit can lead not only to wrong parameter esti- 
mations but also to underestimated parameter errors. To 
handle with this problem, instead of defining a rejection 
parameter we simply performed several fit trials taking 
different wavelength ranges to see how robust a particular 
solution was. In general, there were no significant changes 
in the resulting parameters due to a different choice of the 
fitting interval, thus giving us confidence on the results. 
However, some fits yield \1 < 1. Assuming that the la flux 
errors are not overestimated, such behaviour can occur by 
chance, especially when too few pixels are considered in 
the fit. 

The equivalent-width error estimates are generally 
larger than the propagated error if the flux values are in- 
tegrated along a line. This is because the former errors 
consider the correlation between all three free parame- 
ters. This is not the case for multicomponent fits, where 
the errors in equivalent width may be overestimated since 
the fit algorithm does not consider that the strengths and 
widths of neighbour lines in blends are anti-correlated. 

Concerning the line widths, note that several FWHM 
values in Tables |l| and [2] are slightly larger than the nom- 
inal width of the line spread function (LSF), 2 A. This is 
caused by the difficulty in resolving absorption structures 
separated by less than ~ 200 km s _1 , since single-line fits 
of closely lying line complexes (blends) will inevitably lead 
to large widths. On the other hand, in few cases the mea- 
sured lines are narrower than the LSF. While this might 
partly be caused by a too low placement of the continuum 
near emission lines and BAL troughs in the B spectrum, 
it is certainly not the case for all the narrow lines in the 
A spectrum. In the latter case fitting weak lines also leads 
to low FWHM values (with large associated errors), thus 
reflecting the limitations of the \ 2 minimization using too 
few pixels. 

BAL profiles in the B spectrum have been fitted with 
multicomponent Gaussians to better constrain their equiv- 
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Table 1. Absorption Lines in HS 1216+5032 A. Table 2. Absorption Lines in HS 1216+5032 B. 



Line A oba <y x FWHM W\ a w ID z note 



1 


2250 


.88 


0. 


44 


4 


.99 


1 


.27 





.27 


Hi 1215 





.85155 




2 


2261 


.61 





12 


2 


.10 





.81 





.15 


H I 972 


1 


.32548 




3 


2270 


72 





.19 


4 


.31 


1 


.42 





.20 


Hi 1215 





.86787 




4 


2306 


20 


0. 


16 


2 


.22 





.58 





.19 


H I 1215 





.89706 




5 


2308 


44 


0. 


31 


6 


.66 


3 


.29 





.32 


H I 1215 





.89890 




6 


2331 


.07 





16 


2 


.64 





.75 





.14 


H I 1215 





.91752 




7 


2342 


94 


0. 


17 


4 


.21 


1 


.42 





19 


Fen 2344 


-0 


.00054 


1 


8 


2359 


20 


0. 


15 


2 


.45 





.75 





14 


H I 1215 





.94066 




9 


2363 


08 


0. 


.14 


2 


.41 





.79 





.15 


Hi 1215 





.94385 


10 


10 


2373 


10 


0. 


14 


2 


.46 





.77 





13 


H I 1215 





.95210 




11 


2381 


.84 


0. 


10 


2 


.26 





.99 





13 


Fen 2382 


-0 


.00039 


1 


12 


2385 


30 


0. 


10 


2 


.32 





.99 





14 


Hi 1025 


1. 


.32548 





13 


2392 


.87 


0. 


19 


2 


.61 





.57 





13 


Si iv 1393 





.71690 




14 


2399 


.75 


0. 


45 


5 


.76 





.84 





23 


Hi 1215 





.97402 


2 


15 


2410 


50 


0. 


28 


6 


.18 


2 


.28 





23 


Hi 1215 





.98285 




16 


2412 


.67 


0. 


11 


1 


.58 





.43 





.12 


Hi 1215 





.98464 




17 


2470 


24 


0. 


08 


3 


.38 


1 


.82 





12 


Hi 1215 


1. 


.03200 


3 


18 


2487 


.68 


0. 


13 


1 


.69 





.46 





11 


Hi 1215 


1. 


.04635 




19 


2492 


54 


0. 


21 


2 


.09 





.46 





.15 


Hi 1025 


1. 


.43003 




20 


2509 


24 


0. 


28 


3 


.39 





.55 





15 


H I 1215 


1. 


.06408 


4 


21 


2513 


90 


0. 


14 


2 


.66 


1 


.60 





.25 


Hi 1215 


1. 


.06791 




22 


2536 


43 


0. 


12 


2 


.07 





.49 





09 


Hi 1215 


1 


.08644 


5 


23 


2556 


.38 


0. 


32 


4 


.46 





.86 





19 


Hi 1215 


1. 


.10286 




24 


2561 


73 


0. 


18 


3 


.45 


1 


.04 





.15 


Hi 1215 


1. 


.10726 


7 


25 


2585 


.71 


0. 


11 


2 


.03 





.64 





11 


Fen 2586 


-0. 


.00036 


1 


26 


2591 


15 





12 


2 


.02 





.65 





11 


Hi 1215 


1. 


.13146 




27 


2599 


.35 


0. 


07 


1 


.98 





.88 





10 


Fen 2600 


— 0. 


.00032 


1 


28 


2626 


.59 


0. 


08 


2 


.27 





.89 





.10 


H I 1215 


1. 


.16062 




29 


2647 


.53 


0. 


.17 


1 


.74 





.34 





10 


Cn 1334 





.98387 




30 


2658 


10 


0. 


12 


3 


.17 


1 


.77 





18 


Civ 1548 





.71690 




31 


2662 


.31 


1. 


00 


3 


.01 


1 


.20 





.83 


Civ 1550 





.71677 




32 


2664 


.83 


1. 


00 


2 


.88 


1 


.11 





.81 


Civ 1548 





.72125 




33 


2668 


.68 


0. 


11 


2 


.02 





.84 





15 


Civ 1550 





.72087 




34 


2671 


66 


0. 


20 


2 


.47 





.62 





.15 


Civ 1548 





.72566 




35 


2676 


23 


0. 


23 


1 


.92 





.28 





11 


Civ 1550 


0. 


.72574 




36 


2727 


.61 


0. 


.16 


1 


.80 





.35 





10 


Hi 1215 


1. 


.24371 




37 


2748 


.87 


0. 


22 


4 


.18 





.87 





.15 


Hi 1215 


1. 


.26119 




38 


2791. 


18 


0. 


11 


1 


.28 





.30 





08 


Hi 1215 


1. 


.29600 




39 


2795 


.61 


0. 


07 


3 


.23 


1 


.88 





13 


Mgn 2796 


-0 


.00027 


1 


40 


2801 


.78 


0. 


09 


3 


.51 


2 


.83 





18 


Mgii 2803 


-0. 


.00063 


1 


41 


2805 


50 


0. 


27 


2 


.89 





.73 





.17 


Hi 1215 


1. 


.30778 




42 


2826 


96 


0. 


04 


2 


.43 


1 


.77 





.09 


Hi 1215 


1. 


.32543 




43 


2848 


74 


0. 


15 


1 


.76 





.41 





10 


H I 1215 


1. 


.34335 




44 


2851 


83 


0. 


24 


2 


.52 





.44 





13 


Hi 1215 


1. 


.34589 




45 


2906 


.67 


0. 


11 


2 


.27 





.64 





.10 


H I 1215 


1 


.39100 


9 


40 


2910 


.77 


0. 


2G 


2 


.68 





.35 





11 


Hi 1215 


1. 


.39437 


9 


47 


2927 


54 


0. 


06 


4 


.24 


2 


.33 





10 


Hi 1215 


1. 


.40817 




48 


2939 


28 


2. 


20 


4 


.83 





.64 





52 


Civ 1548 





.89852 


11 


49 


2943 


34 


2. 


05 


4 


.57 





.59 





.53 


Civ 1550 





.89799 


11 


50 


2954 


28 


0. 


05 


2 


.38 


1 


.11 





07 


H I 1215 


1. 


.43016 




51 


2959 


27 


0. 


12 


2 


.08 





.38 





.07 


Hi 1215 


1. 


.43427 




52 


2971 


.79 


0. 


06 


2 


.83 


1 


.04 





.07 


Hi 1215 


1. 


.44457 




53 


2980 


.60 


0. 


18 


1 


.70 





.25 





.07 


Hi 1215 


1. 


.45182 


8 


54 


2982 


93 


0. 


17 


1 


.85 





.34 





11 


Hi 1215 


1. 


.45373 


8 


55 


2985 


43 


0. 


17 


2 


.76 





.84 





.11 


Hi 1215 


1. 


.45579 


8 


56 


3037. 


.03 


0. 


38 


1 


.78 





.11 





.07 


Civ 1548 


0. 


.96166 




57 


3041 


.81 


0. 


53 


1 


.87 





.09 





.08 


Civ 1550 





.96148 




58 


3072 


14 


0. 


15 


1 


.46 





.25 





08 


Civ 1548 





.98434 




59 


3076 


.81 


0. 


23 


1 


.15 





.11 





.07 


Civ 1550 





.98405 




60 


3146 


.54 


0. 


58 


4 


.98 





.49 





18 


Civ 1548 


1. 


.03239 




61 


3151. 


.91 


0. 


19 


1 


.13 





.15 





08 


Civ 1550 


1. 


.03248 




62 


3202 


.78 


0. 


56 


2 


.44 





.17 





.12 


Civ 1548 


1. 


.06872 




63 


3207. 


82 


0. 


70 


2 


.54 





.15 





13 


Civ 1550 


1. 


.06853 




64 


3254 


.73 


0. 


59 


4 


.53 





.47 





20 


Civ 1548 


1. 


.10227 




65 


3260 


.87 


0. 


24 


1 


.84 





.30 





12 


Civ 1550 


1 


.10275 





Line 


A„ b 




"\ 


FWHM 


W x 
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1 Absorption in the ISM. 

2 Or Siiv A1402 at z = 0.721. 

3 Blended with Hi A1025 at z = 1.4082 if b > 30 km s _1 . 

4 Or Sin A1260 at z = 0.984. 

5 Or Sim A1206 at 2 = 1.102. 

6 Blended with Hi A1215 at z = 0.9616. 

7 Or Sin A1260 at z = 1.032. 

8 Associated system. Ly/3 is blended with line A21. 

9 Lines mimic Civ doublet but no corresponding Hi A1215 line is 
detected. 

10 Small contribution from Hi A972 at z = 1.430 to W\. 

11 Strong blend. The identification is supported by the detection of H I 
(A4 and A5) 



1 Absorption in the ISM. 

2 Probably blended with Ly-a. 

3 Blended with Lya at z = 0.984. 

4 Blended with Hi A1025 BAL. 

5 Also S iv A1062 BAL because the line is too strong compared with 
IS Fen A2600. 

6 Or Hi A1025 at z = 1.4082 if b > 30 km s" 1 . 

7 Lines mimic Civ doublet but no corresponding Hi A1215 line is 
detected. 

8 The fit procedure is not able to resolve the second doublet lines; 
therefore, the doublet ratios are probably not real. 

9 Associated system. Ly/3 falls in the Ovi BAL troughs. 

10 Blended with Cn* A1035. Sec § pi for the classification as BAL sys- 
tem. 

11 No line at 3(7 can be identified with Ly/3. 
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alent widths, but these fits do not necessarily represent 
the actual nature of the BAL profiles (see § ||). For con- 
sistency, Tables |l| and || include the Gaussian parameters 
of BAL profiles. 

Absorption lines were identified manually, using vac- 
uum wavelengths and oscillator strengths taken from 
Verner et al. (1994). The identification of lines is not easy 
because of the low resolution and the lack of optical data; 
but it is especially difficult in the B spectrum because nar- 
row intergalactic lines are hidden among the BALs. We 
looked first for interstellar absorption lines in both spec- 
tra. The Mgn A2796,2803 doublet and the Fen lines at 
2586 A and 2600 A are detected in both spectra; the Fe II 
lines at 2344 A and 2382 A are detected at 3cr only in the 
A spectrum (while in the B spectrum Fen A2382 would 
appear at the position of BAL C in A977) . The next step 
was to search for C iv systems and their associated Lya 
and prominent meta l lines. In addition, two Mgn systems 
at z = 0.04 (§ fO|) and 1.14 could be identified in the 
spectrum of B through the positions of the doublet lines 
and their relative strengths. Due to their secure identifi- 
cation, Civ A1548, 1550 and Mgn A2796,2803 lines with 
SNR < 3 have been retained in Tables [j] and & which oth- 
erwise contain only lines with SNR > 3. Lya lines with 
redshifts of 1.16 < z < z e were identified through the cor- 
responding Ly/3 and Ly7 lines, if present. Lines A45, A46, 
B37, and B38 mimic a Civ doublet at z = 0.88, but since 
no significant absorption by Hi is detected at this redshift, 
they have been counted as Lya lines. The remaining lines 
have been considered to be Lya lines and will be discussed 



4. Metal absorption systems in HS 1216+5032 A 
and B 

We now give a description of the most remarkable metal 
absorption systems observed in HS 1216+5032 A and B: 
a Mgn system at z = 0.04, and a strong Civ system at 
z = 0.72. 

4.1. The Civ systems at z = 0.72 in HS 1216+5032 A 
and B 

Three strong C iv systems at z = 0.717, 0.721 and 0.726 
are identified through the AA1548, 1550 doublets at A ~ 
2660 A in both spectra of HS 1216+5032. Although some 
of the C iv lines are under the 3a significance level, the 
well-determined redshifts of the doublet components make 
their identification unambiguous; consequently, all six C IV 
lines appear in Tables @ and g 

The rest frame equivalent widths of the A1548 lines 
range between 0.36 and 1.03 A in the A spectrum, and 
between 0.72 and 1.19 A in B. The large column densities 
implied by these line strengths, make this system a firm 
candidate f or a L yman-limit system (Sargent, Steidel & 
Boksenberg 198E). Also, the large line widths in A and B, 



FWHM > 2 A, suggest that these C IV systems will reveal 
several narrower components at higher resolution. 

No further metal lines are found at these redshifts. Line 
A13 could be identified with Siiv A1393 at z = 0.717, but 
two strong Lya lines at A = 2410 A make the detection 
of the second doublet line impossible. No transitions by 
low-ionization species (e. g., C n A1334) are detected. In 
the low-resolution optical spectrum of HS 1216+5032 B 
(see Fig. [3] in Hagen et al. 1996| ), a significant absorption 
feature at A = 4816 A could tentatively be identified with 
a Mgn A2796,2803 doublet associated to these Civ sys- 
tems, but no absorption feature is seen at this wavelength 
in the optical spectrum of A. 

The small velocity differences of Ava-b ~ —68 ± 28 
km s _1 , —52 ± 31 km s _1 , and —116 ± 30 km s _1 at z = 
0.717, 0.721 and 0.726, respectively, between the C iv lines 
in A and the corresponding ones in B suggest strongly that 
these C iv systems are physically associated. Therefore, 
the large velocity span of roughly 1500 km s _1 along each 
LOS can only be explained if the gas is virialized by a 
cluster of galaxies. The virial mass within the radius given 
by the half separation between the LOSs is derived to be 
~ 10 13 Mq. Consequently, if Civ systems arise in the 
extended halos of galaxies (Bergeron & Boisse 1991), the 
present data gives evidence that the C iv absorbers at z = 
0.72 in HS 1216+5032 A and B arise either in the highly 
ionized intra-group gas of a galaxy cluster, or in small 
structures associated with the cluster. 

The whole absorption complex shows an asymmetry 
in the sense that the weakest system in spectrum A has 
the same redshift as the strongest in B: the ratio of the 
A1548 line's equivalent width in A to that in B varies be- 
tween 1.4 and 3.3. If C iv systems in A and B with similar 
redshifts are physically associated, then the different line 
strengths have direct implications for the size of these ab- 
sorbers. This is because such gas inhomogeneities imply 
that the LOSs sample the clouds on spatial scales similar 
to the transverse cloud sizes. The interpretation remains 
valid if the lines in the present FOS spectra result from 
many unresolved velocity components, because in that 
case stronger lines result from more narrow components 
than weaker lines do, a picture that is still compatible 
with gas inhomogeneities. All these three C IV absorption 
systems should therefore arise in clouds with character- 
istic transverse lengths of ~ 75 ft^kpc, larger than the 
statistical lower limits of ~ 30 /i^kpc derived by Smette 
et al. ( 1995|) for Civ absorbers at z <; 1.5. 

Alternatively, it is still possible that the C iv absorp- 
tion is correlated in redshift, but occurs in distinct, sepa- 



rated structures. At higher resolution, Rauch (1997) has 
shown that in gas associated with metal systems density 
gradients on sub-kpc scales are not uncommon. The logi- 
cal conclusion would be that C IV absorbers are composed 
of a large number of small cloudlets. However, if that is the 
case of the z = 0.72 absorbers, the correlation in redshift 
between lines in A and B is difficult to explain without in- 
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voking cloudlets aligned along a filamentary or sheet-like 
structure. 



4.2. A Mgn System at z = 0.04 in HS 1216+5032 B? 

There is a strong and blended feature at A ~ 2920 A in the 
B spectrum for which we do not find any plausible identi- 
fication with other observed metal systems. One possible 
identification is absorption by two Mg 11 AA2796, 2803 dou- 
blets at z = 0.043 and z = 0.044. Although the Gaussian 
profile fit is not able to resolve the single lines in the red 
trough of the absorption feature, the identification is sup- 
ported by a good match between the line positions and 
profiles of the doublet. The red wing of the whole absorp- 
tion feature could be blended with a Lya line because it 
coincides quite well with a strong Lya line at z — 1.4 in 
A (A47). However, no Fen lines are detected at this red- 
shift, nor further transitions by low-ionization species. The 
Mgi A2852 line is possibly present, but the match in wave- 
length with line B50, AA « 1.4 A, is not very good. Un- 
der the Mg 11-hypothesis, the total rest-frame equivalent 
width of the stronger Mgn doublet component would be 
Wq = 3.97 A, which is somewhat larger than typical values 
found in gas associated with damped Lya (DLA) systems 
at high re dshift (Lu et al. |1996[ ) or in the Milky Way (Sav- 
age et al. 1993| ) . This could be explained if the lines are 
indeed made up of several narrower components, as usu- 
ally seen in DLA systems. The redshift difference between 
both Mg 11 systems implies a velocity span of roughly 300 
km s _1 , also typical of DLA systems. Therefore, there is 
some evidence that these Mgn lines might be associated 
with a DLA system at z — 0.04. 

Regardless of the DLA-system interpretation, however, 
even the absence of Fe 11 and - possibly - also Mg 1 lines 
associated with Mg 11 of such a strength is difficult to ex- 
plain considering the incidence of the former ions in Mg 11- 



selected samples (e.g., Bergeron & Stasinska 1986; Steidel 
& Sargent 1992| ). Consequently, an alternative identifica- 
tion of the A ~ 2920 A feature as Hi A1215 BAL at 
v ~ —6 000 km s _1 (sec §|| and Figure |^) must be con- 
sidered as well. The absence of a corresponding H 1 A1025 
BAL profile, however, is also in this case remarkable, but 
could be explained if the absorber does not completely 
cover the continuum source (see below). 



5. The BAL systems in HS 1216+5032 B: Some 
Qualitative Inferences 

Three BAL systems are observed in the UV spectrum of 
HS 1216+5032 B (see Fig. g, where the system redshifts 
have been arbitrarily numbered 1, 2, and 3). Absorption 
by H 1, C in, N in, and possibly S iv is observed in at least 
two of the systems, while O vi and N v are present in all 
three systems. C11 is observed only in system 2, and Nm 
in systems 2 and 3. 



Many features distinguish these systems from the 
most commonly observed BAL systems (e.g., Turnshck et 
al. 1996j ): (1) the lines are particularly weak and several 
line profiles are not distorted by/or blended with other 
BALs; (2) the maximum outflow velocity v ~ 5 000 km s _1 
is small compared with typical BALQSOs, which exhibit 



terminal velocities of several 10 4 km s _1 (Turnshck 1984) 
(3) C II is present (see Wampler et al. 1995| for other 
BALQSO spectrum with singly-ionized species; Arav et 
al. 1999b ); (4) the strength of Hi in systems 2 and 3 de- 
creas es mo re slowly at the red edge of the troughs (Turn- 
shek |1984|) . 

Also remarkable is the undulating shape of the ab- 
sorbed continuum for A < 2700 A (see Fig. |l|). The po- 
sition of the flux depressions coincides quite well with 
the blue wing of expected emission lines by C ill, N III, 
Ly/3, and O vi. This coincidence seems to suggest that the 
absorption is indeed dominated by very broad line com- 
ponents, which determine the continuum shape, superim- 
posed to the narrower components, here labeled as systems 
1, 2 and 3. The same effect, although less remarkable, is 
observed for S iv and S vi, thus giving evidence for absorp- 
tion by these ions associated with the BAL p henome non 
(also reported for another QSO by Arav et al. 1999a ). 

It is customary to define BALs as a continuous absorp- 
tion with outflow velocities larger than ~ 3 000 km s _1 



from the emission redshift (Weymann et al. 1991) to make 
a distinction between BAL systems and "associated sys- 
tems". However, as pointed out by Arav et al. (1999b), 
such definition does not hold any physical meaning. In 
our case, systems 1 and 2 should be then classified as as- 
sociated systems but, due to our poor resolution, it is not 
possible to establish whether the line profiles are produced 
by continuous absorption or whether they are made of sev- 
eral narrower velocity components (as observed in associ- 
ated systems). The metal lines in these systems (e.g. the 
N v and O vi doublet lines) are relatively narrow and their 
widths could be dominated by the instrumental profile (for 
Hi in system 2, however, the situation is less clear). There- 
fore, the classification of systems 1 and 2 as BAL systems 
must be considered solely instrumental. 

High-resolution spectra of BAL QSOs (Hamann et 



al. 1997) show that BAL profiles do not necessarily result 



from an ensemble of discrete narrow, unresolved lines, but 
as part of a mixture of a continuously accelerated out- 
flow and overlapping narrow components with different 
non-thermal velocity dispersions. For the systems in HS 
1216+5032, it can be assumed that a similar mixture of 
line widths is present, with the broad components domi- 
nating over the narrow ones. In that case, these line pro- 
files should not look so different at higher resolution (with 
exception, maybe, of systems 1 and 2). 

Ionization models have shown that BAL clouds span 
a range of densities a nd/or distances from the ionizing 
source (Turnshek et al. 1996 ; Hamann et al. 1995 ). Fig. [| 
shows that, while the high-ionization species N v and O vi 
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Fig. 2. Broad absorption lines in HS 1216+5032 B plotted 
in velocity relative to z e = 1.451. 

are present in all three systems in HS 1216+5032, Hi, 
C II and doubly-ionized species are absent in system 1 . A 
possible explanation is that the ionization conditions could 
change with outflow velocity, with higher ionization level 
for redshifts closer to z e . However, a quantitative study 
of the ionization conditions in these BAL clouds with the 
present data is made difficult by our inability to reliably 
establish the continuum level at the BAL troughs and thus 
determine column densities (besides the fact that nonblack 
saturation of the line profiles might be present; see Arav 
et al. p-999afl . 

Nevertheless, a more quantitative study of the BAL 
phenomenon in HS 1216+5032 B should be possible us- 
ing higher resolution HST spectra to better estimate the 
QSO continuum and to identify spurious lines at the BAL 
troughs. In particular, BAL system 2 can be studied in 
more detail because it shows more ions and less contami- 
nation by narrow absorption lines at the lines lying in the 



Lya forest. On the other hand, alone medium- resolution 
optical spectra should considerably improve our knowl- 
edge of the BAL systems in HS 1216+5032 through the 
line profiles of C IV, C III, and possibly Mg II, which we sus- 
pect to be present given the presence of C II. The width of 
the trough in system 2 is sufficiently small to study non- 
black saturation effects via doublet ratios provided the 
emission line profiles can be determined (e.g., at the po- 
sition of the Civ BAL). Thus, corrected column densities 
for system 2 appear as very suitable to reliably examine 
photoionization models because the radiation fields can be 
constrained by the presence of low-ionization species. 

6. Lya absorption systems 

We now describe the Lya forest lines observed in the spec- 
tra of HS 1216+5032 A and B. The projected separation 
between LOSs samples well the expected cloud sizes of 
hundreds of kpc if the QSO pair is real (see next subsec- 
tion). However, the effective redshift path where Lya lines 
can in principle be detected is blocked out by the broad 
absorption lines in the B spectrum, so the line samples 
will be small. 

6.1. Lensed or Physical Pair? 

The nature of the double images in HS 1216+5032 is not 
clearly established yet. Maybe the strongest arguments fa- 
voring HS 1216+5032 AB as a physical-pair instead of a 
gravitational lens origin are: (1) BALs are observed only 
in the spectrum of the B image£] (2) no extended, lumi- 
nous object is detected between (or even close to) the QSO 
images down to tur = 22.5; and (3) gravitationally lensed 
QSO pairs normally have image separations 9 ^ 3" (e.g., 
Kochanek, Falco & Munoz 1999| ). In addition, the absence 
of a trend for equivalent widths of intervening Lya lines 
common to both spectra to become similar at redshifts 
closer to z e (see Fig. || and Table ||) also argues for a 
physical-pair; however, this argument is weak as the sam- 
ple of common lines is too small. Unfortunately, for rea- 
sons given in § 2.2, it is not possible to establish whether 



there are differences between A and B in the emission line 
redshifts and line shapes in the spectral region covered by 
the FOS. 

The intriguing point here is that two z ps z e Lya 
systems are detected in both spectra at almost identi- 
cal redshifts to within ~ 90 km s _1 (see § 3.6 for more 
details). Under the physical-pair hypothesis, these ab- 
sorbers should extend over transverse scales greater than 
80 h^kpe. If, on the other hand, HS 1216+5032 AB is a 



2 Curiously, another QSO pair separated by 9", Q1343+266, 
also shows BALs in only one spectrum. On the basis of the 



irregular flux ratio in the optical range, Crotts et al. (1994) 



find evidence that this is a genuine QSO pair . From the optical 
spectra of HS 1216+5032 AB (Hagen et al. |l996| ), however, a 
similar argument is less clear. 
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gravitationally lensed pair, the detection of these systems 
in both spectra can easily be explained by the fact that 
the separation between light paths approaches zero. 



In consequence, the gravitational lens nature of HS 
1216+5032 cannot be yet ruled out. Let us note that in 
a simple lens geometry the virial mass implied by the ve- 
locity span of the C iv absorbers at z 



0.72 (see § (4Jj) is 
consistent with the deflector mass required to produce an 
angular separation between QSO images of 6 = 9". In such 
a model the separation between light paths at the source 
position is derived to be 9 S = 17". Therefore, confirma- 
tion of HS 1216+5032 as a gravitationally lensed QSO 
would have dramatic consequences for our understanding 
of the BAL phenomenon. As pointed out by Hagen et al. 
( |1996|) , if HS 1216+5032 AB were indeed the mirror im- 
ages of one unique QSO, the coverage fraction of BAL 
clouds would have to be very small. Incomplete coverage 
of the continuum source or the BLR in BAL clouds have 
been confirmed using high-resolution spectra (Hamann et 
al. 1997 ; Barlow & Sargent 1997 ). In fact, the cloud sizes 
derived from BAL variability could be 10 4 times smaller 
than the distances to the source derived from photoion- 



ization models (Hamann et al. 1995), implying subtended 
angles not much larger than 9 S . Therefore, the scenario in 
which only one light path crosses the BAL clouds is not 
unrealistic, and the non-BAL spectrum of HS 1216+5032 
A can also be considered if the QSO pair is lensed.^ 



In summary, there are reasons to believe that HS 
1216+5032 is a binary QSO, but a gravitational lens ori- 
gin of the double images cannot yet be excluded. A defini- 
tive assessment of this issue must await medium resolution 
optical spectra to better derive emission redshifts and to 
establish differences in the continuum or emission lines. 
Alternatively, deep infrared or radio observations would 
help determining the nature of the double QSO images. 
The non-detection of the lensing galaxy/cluster would be 
a strong proof for the binary QSO hypothesis, since wide 
separation (9 > 3") "dark lenses" enter in conflict with 
current models of structure formation (Kochanek, Falco 
& Munoz 1999 ), which predict few w ide separation gravi- 
tational lenses (e.g., Kochanek 1995 ). 



Throughout this section we will assume that HS 
1216+5032 AB is a physical pair. The proper separation 
S(z) between LOSs obeys then the relation S = 9 D oc , 
where 9 is the angular separation between the A and B 
QSO images, and D oc is the angular-diameter distance 
between the observer and a cloud at redshift z. For the 
redshift range of interest S ~ 80 h^ kpc. 




0.8 



Fig. 3. Three a rest-frame detection limits of Lya lines 
in HS 1216+5032A (thick line) and B. 

6.2. Line List 

Table || lists Lya lines detected at the 3cr level in HS 
1216+5032 A and B. Figure ^ shows 3cr detection thresh- 
olds as a function of redshift. Out of 36 lines in the spec- 
trum of A, a total of 23 lines within z = 0.85 to 1.43 
have Wo (A) > 0.24 A. This redshift interval excludes 
lines within 3 000 km s _1 from the QSO redshift to avoid 
the "proximity effect" (Bechtold 1994). The derived num- 
ber density is log{dN/dz) = 1.58^°^ at log(l + z) = 
0.33 ± 0.06, in very good agreement with the results of 

The scarce 



1998) 



the HST Key Project (Weymann et al. 
sample of lines in B is clearly explained by the shorter 
redshift path allowed by the BAL profiles. 

6.3. Definition of Line Samples 

To define the number of coincidences and anti- 
coincidences we have selected from all 36 Lya lines ob- 
served in A (the spectrum with better signal-to-noise) 
those ones (1) at z < 1.43; (2) not associated with metal 
lines; (3) at wavelengths not covered by the BAL troughs 
in B; and (4) at wavelengths where 3er detection limits in 
B are lower than the measured equivalent width of the line 
in A. Lines in B that had SNR > 3 and fulfilled criteria (1) 
to (3) were also selected. In some cases, lines in A have a 
corresponding absorption feature at the same wavelength 
in B, but at a too low significance level to unambiguously 
designate the line pair as a coincidence. These lines in A 
were excluded (A26, A36, A44). Notice that criterion (3) 
automatically prevents comparing line equivalent widths 
for which the uncertainties introduced by the placement 
of the B continuum are large. 

The final sample of Lya lines suitable for this study is 
composed of N = 11 redshift systems within 1.16 < z < 
1.43, a range not much smaller than the one allowed by 
conditions (1) to (4). Out of this number, Nq = 6 systems 
show Lya lines in both spectra, and N& — 5 lines are 
detected in only one spectrum. The rest-frame equivalent 
widths range from 0.14 to 0.76 A in A, and from 0.23 to 
1.16 A in B. This sample is hereafter called "full sample". 



3 Let us mention, however, that in the lensed-pair hypothe- 
sis the fact that BAL absorption is produced by many clouds 



(Hamann et al. 1997 ) lessens considerably the probability that 



LOS A does not encounter neither such clouds. 
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Table 3. Lya Lines with W > 3<r w in HS 1216+5032 A and B. 
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1 Three sigma detection limits in B > W(A). 

2 Spectral region covered by the BAL troughs in B. 

3 C IV associated. 

4 Maybe S iv A944 BAL. 

5 Absorption feature present in B at low significance. 

6 Associated system. 



Coincident lines are shown in Figure [|, where the thick 
line represents the flux of A and the zero velocity point 
corresponds to the redshift of the A line. Only the first 
six panels show lines in the "full sample"; the remaining 
lines arise in Lya clouds likely to be influenced by the 
QSO flux. All coincident lines are separated by less than 
100 km s _1 . Anti-coincident lines are all uniquely defined, 
as there are no corresponding lines in the other spectrum 
within several hundred km s . 

However, there are two exceptions that must be 
pointed out. The first one is line 29 in B, separated by 
198 km s -1 from the nearest Lya line in A (37). These 
lines in A and B have been counted as anti-coincidences, 
though we are suspicious of this interpretation as the line 
profiles suggest line B29 is blended with an absorption 
feature seen in both spectra. The second is another of 



the anti-coincidences, line B50, which could alternatively 
be identified with Mgi at z = 0.04. These cases will in- 
troduce unavoidable uncertainties in the results presented 
here. Let us recall that the exclusion of any one line from 
the samples would lead to significantly different results 
for R c . This illustrates the real uncertainties dominating 
simulations with such a small number of observed lines. 
Moreover, the samples are limited by the large redshift 
path blocked out by the BAL troughs in B, so there must 
be a considerable loss of information. 

A line significance level of 5 for lines in A has been cho- 
sen to perform the likelihood analysis described in the next 
subsection . This selection automatically excludes two of 
the redshift systems in the full sample and defines the 
following sub-samples: 
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Fig. 4. Lya lines present both in HS 1216+5032 A and 
B, with SNR > 3, Av(A - B) < 100 km s -1 , and no met- 
als associated. These lines lie at wavelengths not covered 
by the BAL troughs in B. The right number indicates to 
which sample the line pair belongs. The three bottom pan- 
els show lines that fall within 3 000 km s _1 from the QSO 
redshift. The thick line represents the flux of A. The ab- 
scissa is in km s _1 and the zero- velocity point is defined 
by the redshift of each A line. 



2. Strong+weak line sample (sample 2): same as previous 
sample but with Wq > 0.17 A. This selection yields 
N c = 4 and N A = 5. 

6.4- Maximum Likelihood Analysis 

In the following we discuss a likelihood analysis on cloud 
sizes in front of HS 1216+5032. The technique is based 
on the definition of a likelihood function C that gives the 
probability of getting the observed number of coincidences 
and anti-coincidences. Evidently, such a function must de- 
pend on the shape of the absorber. Here we concentrate on 
two possible cloud geometries for which C can be derived 
analytically: spheres and cylinders. 

For coherent spherical clouds, the probability that one 
cloud at redshift z is intersected by one LOS is given by 
McGill QL990D as: 



<j>(X) = - jarccos [X{z)\ - X(z)y/l - AT(z) 2 | , (1) 

where X(z) = S(z)/2R C and R c is the cloud radius. The 
samples consider only redshifts where at least one LOS in- 
tersects a cloud (otherwise one should make assumptions 
on the cloud distribution along the LOS). Therefore, one 
must compute the probability that, given that a line ap- 
pears in one spectrum, a line will appear in the other spec- 
trum. This probability is given by Dinshaw et al 
as: 



(1997) 
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Fig. 5. Lya lines present either in HS 1216+5032 A or B, 
with SNR > 3, Aw (A - B) < 100 km s _1 , and no met- 
als associated. These lines lie at wavelengths not covered 
by the BAL troughs in B. The right number indicates 
to which sample the line pair belongs. The abscissa is in 
km s _1 and the zero-velocity point is defined by the red- 
shift of each present line. 



1. Strong-line sample (sample 1): lines for which SNR(A) 
> 5, SNR(B) > 3, and W > 0.32 A. This selec- 
tion yields Nq = 2 coincidences and N& — 4 anti- 
coincidences. 



Finally, the probability of getting the observed number of 
coincidences and anti-coincidences is given by the product 



(3) 



where the indexes i and j number the coincidences and 
anti-coincidences, respectively. Note that if there is at least 
one anti-coincidence, then C has a maximum value. Oth- 
erwise it grows monotonically. 

The solid curve in Fig. ^ shows the results of the like- 
lihood function C(R C ) normalized to its peak intensity for 
the various samples. The figure also shows the cumula- 
tive distribution of C(R C ). The estimated 2a limits on 
cloud diameters for HS 1216+5032 - derived from the cu- 
mulative distribution - are listed in Table | The peak 
intensity of C is used to find the most probably radii, 
i.e., 96 h^Q kpc for the strong line sample and 128 
kpc for the strong+weak line sample. From the cumula- 
tive distribution, the respective median values are 145 and 
173 kpc. In both cases, the derived sizes are larger if 
weaker lines are used. This result provides evidence that 
Lya clouds must have a smooth density distribution. 

In a second model, let us suppose that Lya lines occur 
in filamentary structures lying perpendicular to the LOSs. 
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Fig. 6. Likelihood function C(R C ) normalized to its peak 
intensity vs. cloud radius R c , and cumulative distribution 
for spherical (solid line) and cylindric absorbers. 

Such structures can be idealized as cylinders with a radius- 
to-length ratio <C 1. For cylindric clouds, the probability 
that one LOS intersects a cloud at redshift z given that 
the other LOS already does, reads (A. Smette, private 
communication) : 

f jarcsin [X(z)- 1 ] 

+ ^/X{zf - 1 - X(z)\ forX(z) > 1; (4) 



1 - ^ X(z) 



otherwise, 



where X(z) = S(z)/2R C and R c is the radius of the cylin- 
der. The dotted curve in Fig. ^ shows the results of the 
likelihood function C(R C ) normalized to its peak inten- 
sity and the cumulative distribution for the various sam- 
ples. Most probably cylinder radii are 49, 65, and 85 h§Q 
kpc for the strong, strong+weak and full line samples, re- 
spectively. The respective median values are 70, 87 and 
115 kpc. Two sigma bounds are displayed in Ta- 

ble ^. Again, the derived sizes are larger if weaker lines are 
used, providing evidence that the density distribution in 
Lya clouds must be smooth (see, e.g., Monier et al. 1999| ). 
These size estimates are almost 50% lower than for spher- 
ical absorbers. This can be explained by the fact that, in 
general, elongated structures can more easily reproduce 
the observed number of coincidences for a given radius 
than a sphere; in fact, the probability 4>(R C ) for spheri- 
cal absorbers vanishes for R c < S/2, while for cylinders 
4> > for all nonzero values of R c . The lengths of such 
structures are much larger than these values, but unde- 
fined in the model. 



6.5. Equivalent Widths 

The sample of Lya lines common to both spectra gives 
model-independent information about the size scales of 
the absorbers. Figure ^ shows the rest-frame equivalent 
widths of lines common to both spectra and their la er- 
rors, where the dashed straight line has a slope of unity. 
Only lines without associated metal lines are included. 
Note that this selection yields line-pairs with A > 2620 A, 
so the influence of uncertainties in the continuum place- 



Table 4. Diameter D of spherical and cylindric Lya 
clouds in HS 1216+5032 in the redshift interval z = 
1.15 — 1.43 as derived from the likelihood function defined 
in (3). Values are in kpc and represent 2a limits. Hq = 50 
km s _1 Mpc" 1 and go = 0.5. 



Cloud shape 


Samplel 


Sample2 


Sample3 


Spherical 


136 < D < 880 


172 < D < 896 


224 < D < 1160 


Cylindric 


52 < D < 454 


86 < D < 460 


112 < D < 610 



Sample 1: Lines with Wo > 0.32 A; 2 coincidences, 4 anti- 
coincidences. 

Sample 2: Lines with Wo > 0.17 A; 4 coincidences, 5 anti- 
coincidences. 

Sample 3: Full sample: 6 coincidences, 5 anti-coincidences. 



mcnt (e.g. due to the BAL troughs or to spectral regions 
with high noise level) over the equivalent-width estimates 
is minimized. We have included one line pair within 3 000 
km s- 1 from the QSO's redshift (A52,B49) but not the 
"associated system". The dotted crosses represent line 
pairs for which the B line has 2.58 < SNR < 3, i.e., those 
ones labeled with footnote 5 in Table ||. Additionally, to 
illustrate the significance of the non-detections we have 
also included lines detected in only one spectrum, plotted 
against the corresponding 3a detection limit in the other 
spectrum. 

Most of the coincident lines are stronger in the B spec- 
trum, which is explained by the difficulty in detecting lines 
in B (an additional consequence of this is that most anti- 
coincidences are lines in B). It can be seen that there are 
some significant deviations from AWq = 0, even exclud- 
ing the line pair within 3 000 km s . For the latter such 
deviation would be expected, since the QSOs themselves 
might be an important ionizing agent in clouds next to 
the QSOs.0 

The equivalent width differences AWq — \Wo(A) — 
Wo(B)| show no correlation with Av. Instead, AWo and 
max(FFo(A), Wo(B)) seem to be correlated (see Fig. ||), 
with larger equivalent width differences for larger equiv- 
(e. 



alent widths (e.g. Fang et al. 1996 ). Both effects, al- 
though not statistically significant, suggest coherent struc- 
tur es, i.e ., no "cloudlets" at similar redshifts (Charlton et 
al. 1995 ). In conclusion, LOS A and B must sample co- 
herent clouds showing small - but otherwise significant - 
density gradients on spatial scales of ~ 80 /i^kpc, the 
linear separation between LOSs in this redshift range. 



4 Indeed, if the QSO radiation field contributed significantly 
to ionize these clouds, one would expect 14^ (A) < W^B) for 
these lines, because QSO A is more luminous than B. Such 
difference in equivalent width is, however, not observed. 
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Fig. 7. Rest-frame equivalent widths of Lya lines com- 
mon to both spectra, with SNR > 3, Av < 100 km s _1 , 
and no metal lines associated. Dotted crosses represent 
line pairs with 2.58 < SNR < 3. The "associated system" 
has not been considered. Lines detected in only one spec- 
trum are also shown, with upper limits representing 3er 
detection limits in the other spectrum. 



6.6. The z a w z e Lya systems in HS 1216+5032 A and B 

Two z a w z e Lya systems are observed in the spectra 
of both HS 1216+5032 A and B, through the partially 
resolved Lya lines A53, A54, A55, and B51 and B52 (see 
Fig. ||). The lines in A and B differ by Aua53-B51 = -51 
km s _1 and Ai>A54-B52 = —87 km s . The system in A 
is blueshifted by 122 km s~ x relative to z e (A), while the 
system in B is redshifted by 280 km s^" 1 relative to z e (B) 
(but see § 2.2). Line A55 has no clear counterpart in B. 



No metal lines are found associated with these systems. 
The total equivalent width of the A lines is larger than in 
B by a factor of 2.1. 

The absence of highly ionized species such as N V and 
Ovi suggests that these systems are probably not ph ysi- 
cally associated with the QSOs (e.g., Turnshek 1984; Pe- 
titjean 1994). Furthermore, the small velocity differences 
between lines in A and B, and the fact that both sys- 
tems have two line components well correlated in redshift, 
strongly suggest they arise in the same absorber (see also 



Petitjean et al. 1998; Shaver & Robertson 1983). If this is 



true, the different Lya line strengths between A and B im- 
ply either characteristic transverse sizes of - 80 h^kpc 
or maybe even larger clouds with ionization conditions 
that change on such scales. The size scales are compatible 
with the idea that these systems arise, for example, in the 



Fig. 8. Absolute value of rest-frame equivalent width dif- 
ferences (A) vs. max(M / o(A), Wo(B)) of Lya lines with 
SNR > 3, Av < 100 km s _1 , and no metal lines asso- 
ciated. The "associated system" has not been considered. 



very extended halo of an intervening galaxy (Weymann 
et al. 1979), or in the intra-group gas of the QSO host 
galaxies. 

If the radiation field of the QSOs is the main ionization 
source in the clouds giving rise to these systems (assuming 
the QSO pair is real), then the systems in B should be less 
ionized than the A ones. This is because QSO A is intrin- 
sically more luminous than B in the ionizing continuum. 
Consequently, the different line strengths in A and B can 
also be explained if the LOSs cross regions of similar den- 
sity, with LOS B probing regions with more neutral gas. 
Of course, the latter case does not rule out the intervening 
nature of these systems. 



6.7. Discussion 



Fang et al. ( 1996) ) have pointed out that there seems to be 
a trend of larger estimated cloud sizes with increasing LOS 
separation S (see their Fig. 5). If such a trend were real, it 
would imply that the scenario of uniform-sized spherical 
clouds is too idealized. Their data, however, are insuffi- 
cient to discern whether the effect is due to non-uniform 
cloud size or simply non-spherical cloud geometry. The re- 
sults presented in our work (slightly modified by defining 
a sample of lines with Wo > 0.40 A to be consistent with 
these authors) fit this tendency well because they add a 
new measurement of relatively small cloud sizes at small 
LOS separation. We note, however, that D'Odorico et al. 
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( |1998| ), who use a larger database and an improved statis- 
tical approach, find no correlation between cloud size and 
LOS separation. Clearly, more QSO pairs are needed that 
span a range in LOS separations to confirm or discard the 
size/separation correlation. 

Here we propose that the notion of a single population 
of uniform-sized clouds must be revised. In fact, if Lya 
clouds span a range of sizes between, say, a few hundred 
kpc to a few Mpc, then LOSs to QSO pairs with arcminute 
separations would be crossing not only huge and coherent 
structures, but also smaller clouds correlated in redshift. 
Further evidence for a more complicated scenario than 
usually assumed comes from cosmological hydrodynamic 
simulations made for z = 3, which show entities of non- 
uniform sizes grouped along filamentary structures (e.g., 
Cen & Simcoe 1997). However, we recall that the situation 
at lower redshift can be different if such structures do 
evolve in size. 

Size evolution of Lya absorbers has not yet been ob- 
served, in part because the number of adjacent QSOs with 
suitable angular separations is not statistically significant, 
but also due to the scarce number of observations at low 
redshift. In particular, the present data on HS 1216+5032 
(< z > — 1.3) do not c onfirm the suggestion (Fang et 
al. 11996; Dinshaw et al. 199S) that cloud sizes increase 



with decreasing redshift. This result is in variance with 



the findings by D'Odorico et al. (1998) 



Assuming that the Lya absorbers are filamentary 
structures - modelled as cylinders of infinite length - ly- 
ing perpendicular to the LOSs, our likelihood analysis 
leads to almost 50% smaller transverse dimensions than 
for spherical clouds. Unfortunately, the method presented 
here is not capable of distinguishing between cylinders 
and spheres. Moreover, the information provided by two 
LOSs might be insufficient to make such a distinction 
possible, so o ne should use triply imaged QSOs (Crotts 
& Fang |1998| ) or, even better, QSO groups (Monier et 
al. |l99S| ). It is worth saying, however, that flattened struc- 
tures are more able to simultaneously reproduce the re- 
quirements of neutral gas density from photoionization 
models and the transverse scale lengths derived from dou- 
ble LOSs than spherical absorbers do. Such photoioniza- 
tion models lead to structures with thickness-to-length ra- 
tios of ~ 1/30 (Rauch & Haehnelt 1995 ). Indeed, hydro- 
dynamical simulations in the context of hierarchical struc- 
ture formation have shown that at z ~ 3, high density gas 
regions (producing Wo > 0.3 A Lya absorption lines) are 
connected by filamentary and sheet-like structures roughly 
10 6 times less dense than the embedded condensations. 
The filaments seem to evolve slowly and still fill the Uni- 
verse at z <~ 1 (Dave et al. 1999; Riediger, Petitjean & 



Mucket 1998), giving rise to the majority of strong Lya 
lines. If this picture is correct, the requirement of differ- 
ent cloud populations might not be necessary to explain 
current observations. 



New high-resolution ultraviolet observations of QSO 
pairs are needed. Despite technical difficulties (close QSO 
pairs tend to have such different observed fluxes that 
high dispersion spectroscopy normally leads to at least 
one spectrum with high noise levels) they should im- 
prove considerably our knowledge of the Lya-cloud geom- 
etry by (1) analyzing absorbers that produce only weak 
(Wo ~ 0.01 A) lines, i.e., those surely not associated with 
low-brightness galaxies, and (2) testing models that con- 
sider column density distributions. 

7. Summary 

We have presented HST FOS spectra of the QSO pair HS 
1216+5032 AB. Our results are summarized as follows: 

1. Metal systems. Three strong and complex C iv absorp- 
tion systems are observed in both HST spectra at z w 
0.72. The velocity span along the LOSs, Av ~ 1500 
km s _1 , is large and suggests the gas where this Civ 
occurs might be associated with a cluster of galaxies. 
Differences in the line strengths suggest that, if these 
systems arise in coherent structures, they must have 
characteristic transverse sizes of ~ 75 h^kpe. Alter- 
natively, the systems may be composed of a large num- 
ber of cloudlets correlated in redshift. 

In the spectrum of B, a strong absorption feature at 
A ~ 2910 A is identified with two Mgn doublets at 
z = 0.04. If this identification is correct, these lines 
could arise in a low-redshift damped Lya system. 

2. The BAL systems in B: The spectrum of QSO image 
B shows BAL troughs by Hi, Cn, Cm, Nin, Nv, 
Ovi, and possibly Siv and Svi arising in absorption 
systems at outflow velocities from the QSO of up to 
5000 km s" 1 . The BAL troughs arise probably as a 
consequence of absorption by a mixture of broad and 
narrow components. 

3. Lya absorbers. Due to the redshift path blocked out by 
the BAL troughs in B the number of detected lines is 
small. Selecting lines not associated with metal lines, 
with Av(A - B) < 100 km s" 1 and W Q > 0.17 A yields 
four lines common to both spectra and five lines with- 
out counterpart in the other spectrum. For Wq > 0.32 
A lines the numbers are two and four, respectively. 
Using a maximum likelihood technique, most proba- 
bly diameters for spherical clouds of 192 and 256 h^Q 
kpc are found for W > 0.32 A and W > 0.17 A 
lines, respectively. The 2er limits derived using the cu- 
mulative distribution of the probability function are 
136 < D < 880 kpc and 172 < D < 896 

kpc for the respective samples at < z >— 1.3. Assum- 
ing that the absorbers are filamentary structures lying- 
perpendicular to the LOSs, transverse dimensions al- 
most 50% smaller than for spherical clouds are found. 
In both cases, the results of our analysis do not con- 
firm the claim that the characteristic size of the Lya 
absorbers increases with decreasing redshift. 
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Independently of the cloud models used, we note that 
there are significant equivalent width differences be- 
tween lines in A and B. Also, there appears to be a 
trend of larger equivalent width differences with in- 
creasing line strength, while no velocity differences be- 
tween common lines is found. This provides evidence 
that the absorbers are coherent entities. The results for 
each line sample suggest that the absorbers must have 
a smooth gas density distribution, with lower density 
gas being more extended. 
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Fig. 1. HSTFOS G270H spectra of HS 1216+5032 A and B and their la errors. The dotted line represents the fitted 
continuum (§ 2J3) and the ticks mark the positions of 3<r absorption lines (§ ||), with exception of less significant ClV 
and Mgii doublet lines whose identification is supported by the detection of the corresponding Lya line (see § || and 
Tables p] and ||). The vertical bars in the B plots indicate the expected position of emission lines at z e — 1.451. The 
horizontal bars correspond to a velocity of 5 000 km s _1 in the QSO rest frame, a region thought to be affected by the 
BALs (see also Fig. ||). Notice in B the broad troughs on the blue side of the Siv, Ovi, Nm, Cm and possibly Svi 
resonance lines. 



